Depleted pyrochlore antiferromagnets 

Christopher L. Henley 

Dept. of Physics, Cornell Univ., Ithaca NY 14853-2501 USA 
E-mail: clh@ccmr.cornell.edu 

Abstract. I consider the class of "depleted pyrochlore" lattices of corner-sharing triangles, 
made by removing spins from a pyrochlore lattice such that every tetrahedron loses exactly one. 
Previously known examples are the "hyperkagome" and "kagome staircase". I give criteria 
in terms of loops for whether a given depleted lattice can order analogous to the kagome 
"-^/S X \/3" state, and also show how the pseudo-dipolar correlations (due to local constraints) 
generalize to even the random depleted case. 



Consider "bisimplex" antiferromagnets ^\, meaning that every spin is shared between exactly 
two triangles or tetrahedraLet there be only the nearest-neighbour coupling, so the system is 
highly frustrated. After the kagome [2], pyrochlore [3lH|, and garnet f5j lattices, further such 
systems were discovered such as the "kagome staircase" (realised in e.g. Ni3V208 [Sj) and the 
"hyperkagome" lattice [U [7j (realised in Na4lr308 [9]). Both magnetic lattices are obtained 
from the pyrochlore lattice by removing 1/4 of the sites so as to leave a network of corner- 
sharing triangles. This paper considers the entire family of such "depleted" structures, and the 
equilibrium states of classical or semiclassical spins on them. 

The pyrochlore lattice is most simply visualised as the "medial graph" (bond midpoints) of a 
diamond lattice. A depleted lattice is made by placing a dimer covering on the diamond lattice 
bonds, then removing spins from the covered sites. This constraint is plausible: in real spinels 
with magnetic B sites (= a "pyrochlore" magnetic lattice), dilution is achieved by substituting a 
nonmagnetic species X on B sites. Due to size (or maybe charge) imbalance, X ions repel; hence 
the lattice gas of X ions maps to an Ising antiferromagnet [3] in a field. That is itself a highly 
frustrated problem, with a degenerate ground state (all those dimer coverings) - provided the 
lattice gas has only nearest neighbour interactions. Further neighbour ion terms presumably 
select a specific depletion pattern. 

After some examples of depleted lattices, I address two questions (for the periodic and random 
cases): (i) What is the pattern of magnetic order (if any); (iii) How do we generalize the 
disordered classical liquid with pseudo-dipolar correlations due to the constraints? 



1. Regularly depleted pyrochlore lattices 

In this section, I catalog some highly symmetrical dimer coverings of the diamond lattice, each 
of which specifies a different depleted pyrochlore lattice. 

^ I will not consider other depleted lattices, except to observe the hidden symmetry (considering nearest- 
neighbour bonds) whereby the hyperkagome lattice ^ is equivalent to a garnet lattice [5], is reminiscent of 
the hidden equivalence of the 1/5-depleted square lattice of CaV409 _K to the 4-8 lattice [11| . 



These can be conveniently be visualized in two possible ways: (i) projecting the conventional 
cubic cell (containing four diamond sites, i.e. four pyrochlore tetrahedra) in the (100) direction; 
or, (ii) expressing the diamond lattice as a stacking of puckered honeycomb layers, in which the 
(odd) sites have an additional bond extending upwards (downwards) to the next layer. 

1.1. Cubic conventional cell 

Consider the family of patterns preserving the periodicity of the conventional cubic cell, with 
four dimers per cell. There are only three symmetry-inequivalent ways to place them. Pattern 
1 has all dimers oriented the same; the diamond lattice separates into disconnected (puckered) 
honeycomb layers, and the spins form stacks of uncoupled kagome lattices, as induced by a field 
in "kagome ice" [12j . 

Pattern 2 has two dimers in one orientation and two in another orientation. This again 
separates the lattice into disjoint slabs, now transverse to a (110) axis. The depleted lattice is 
a "kagome staircase", which has the topology of a kagome lattice, but folded so the hexagons 
alternate between two (111) type orientations. Imagine a canonical spinel AB2O4 with cubic 
lattice constant ac ~ 8.2 A; depletion by B-site vacancies gives the formula A1B1.5O4 with 
X = 1/4. Pattern 2 makes the structure orthorhombic, with a ~ ad sqrtl^ b « ^/2ac, and c « Oc. 
This nearly describes the "kagome staircase" compounds e.g. VC01.5O4, except each successive 
slab is slid a half cell relative to the one before. 

Pattern 3 uses each of the four possible orientations once, yielding the "hyperkagome" 
arrangement with cubic symmetry. as realized [9] in the spinel Na4lr308, i.e. 
"(Nai.5)i(Ir3/4Nai/4)204" in our framework, 

1.2. Kagome layer stacking approach 

I will focus on the layered patterns made by viewing the diamond lattice as a stack of puckered 
honeycomb layers (a kagome layer of spin sites) connected by vertical linking bonds (forming a 
triangular lattice). Say a fraction xunk of linking bonds is depleted; then the honeycomb layer 
has depletion rc^ag = (1 — xiink)/3; its depletion pattern is a dimer covering having monomers 
at the endpoints of the depleted linking bonds. If either kind of layer lacks threefold point 
symmetry, we restore it by stacking the layers with a rotation (producing a screw axis). 



Table 1. Some depleted lattices, with properties of the shortest loops (Loop density is per site; 
in the tags, interlayer bonds are underlined for comparison to Fig. [H) 
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Fig. [T] shows some examples; there are many more, e.g. in Fig. HJc) a different honeycomb 
dimer pattern could be used. Fig.[T][^b) is just the cubic "hyperkagome" structure. With xn^]^ = 
and Xkag = 1/3 we get a dimer covering of the honeycomb lattice. 




Figure 1. Depleted lattices built by stacking puckered layers along a 3- fold axis. Diamond 
lattice are shown bonds as lines (in the layer) or as circles (for linking bonds: up=shaded, 
down=black). Removed bonds are indicated by dashed lines or empty circles. In (a) and (c), 
successive layers are stacked with a 27r/3 rotation around the dotted triangle, (a), a -v/3 x -v/3 
pattern; (b). hyperkagome lattice; (c). another 2x2 pattern; (d). a \/7 x \fl pattern. 

2. Magnetic ground state 

The magnetic Hamiltonian is assumed to be W = J^i^ij^ Sj • Sj where only nearest-neighbour 
isotropic interactions are included, and {sj} are either classical unit vectors or quantum 
spins with 5 S> 1. Let us review the known story for the kagome or garnet/hyperkagome) 
antiferromagnets. The classical ground states are the (many) configurations in which the spins 
differ by angles 27r/3 on every triangle. It then transpires that (i) a (still highly degenerate) 
subset of "coplanar" states gets selected at harmonic order, which amounts to a 3-colouring; 
(ii) a specific coplanar ordered state is selected by anharmonic fluctuations, which was the 
"\/3 X \/3" state in the kagome case [HI [131 [15] ^-iid Lawler's state in the hyperkagome/garnet 
case [7]. Do these results generalise? 

2.1. Coplanar states 

At harmonic order, thermal fluctuations (in the classical case t2ii4j) or quantum fluctuations [13^ 
[TB] select "coplanar" ground states, such that spins in different triangles lie in the same plane 
of spin space. [□ This selection should carry over to arbitrary corner-sharing triangle networks. 

A coplanar ground state has every spin in one of three directions so it is effectively a 3- 
colouring of the sites (equivalently, diamond-lattice bonds) by colours A, B, C, such that every 
triangle (i.e. bonds meeting at one diamond vertex) has one of each colour. By Konig's 
theorem [T7j of graph theory, a 3-colouring exists on any bipartite graph with coordination 
z = 3, and hence on any depleted lattice. 

Any 3-colouring of a depleted lattice can usefully be reimagined as a 4-colouring of the 
original pyrochlore lattice, with the fourth color corresponding to the depleted sites. Thus, any 
3-colouring in fact generates four different ways to build a depleted lattice (along with a sample 
3-colouring of each), depending on which color we select for removal. (In Lawler's state [TJ, all 
four colours are equivalent; each colour forms a "trillium" lattice [18j.) 

2.2. Ising mapping and effective Hamiltonian 

As is well known, at harmonic order all the coplanar states have equivalent Hamiltonians [2lll3j. 
Hence, fluctuations distinguish among them only at anharmonic order. Observe too that 
coplanar configurations cannot distinguished on the loop-free "cactus lattice" - the medial graph 
of a z = 3 Bethe lattice - since they are all symmetry-equivalent by permutations of the sites. 
Loops are essential to state selection [I9l |20] . 

On the kagome lattice, coplanar states are more transparently represented by "chiralities" 
r]a = ±1, Ising variables defined on triangle centers a, and equal to -|-1 (—1); if the colours 

The constraint counting argument of 4 carries over independent of how the triangles are arranged. 



ABC run counterclockwise (clockwise) around the triangle. The colouring can be uniquely 
reconstructured (modulo trivial symmetries) from {rja}, but not every configuration of {ija} 
corresponds to a colouring. With approximations, one can obtain an effective Hamiltonian of 
form TicS = — Y^aji •^ajSVaVp ill the quantum case [E], and also in the classical case [2T] at 
small T, based on [22]. The nearest neighbor Ising coupling J\ < 0, so the optimum state is an 
antiferromagnetic pattern of ry^ on the honeycomb vertices. 

What about d = 3? Let the index a label triangles, or equivalently diamond-lattice vertices. 
A gauge choice is necessary on every triangle to define which sense of its normal vector is "up" , 
before the spin chirality can be defined. I conjecture that, in general, the sign of Ji is the 
opposite of the projection of the normal vectors of the respective triangles. 

Then the ground state has an alternating chirality pattern (which is always possible, since 
the diamond vertices are bipartite). For the hyperkagome lattice - the most regular depleted 
example - this gives the correct answer: Lawler's state, favored in a large-n calculation [7j and 
found in simulations [8].) Whereas the "-y/S x ^/3" state of the kagome had alternating colours 
(e.g. ABABAB) around loops and triple colors {ABC ABC.) along lines, and has a nonzero 
ordering vector, Lawler's state has alternating colours along lines and ordering vector Q = 0. H 

The basis for believing < is general is that in the kagome case, it was expressed [151 EI] 
in terms of expectations of spin-wave fluctuations of "soft" modes, which have only anharmonic- 
order restoring forces. Such soft modes are generic to the coplanar state on any of our lattices. 
They sum to zero on every triangle, a "zero-divergence" constraint that implies generic power- 
law correlations of the fluctuations, whether classical [23l|2ll[25] or quantum [15l|26]. We obtain 
J7i < if the sign depends on orientation the same way it does asympotically. A caveat is that 
at root, the crucial spin-wave correlations must depend on the loops (as noted above, only loops 
distinguish among coplanar states); the pseudo-dipolar correlations are just a coarse-grained way 
to incorporate the net effects of many long loops. An alternate local derivation of TleS, based 
on a loop expansion [191 EOl EBj might better capture the differences among depleted lattices. 

2. 3. Colourability 

Let's call the lattice "colourable" if there is a 3-colouring satisfying the above rule. That is true 
if and only if every fundamental loop of diamond-lattice bonds is colourable. Tag the each bond 
of the loop as "s" or "6" depending whether the corresponding site is analogous to a "straight" 
or "bent" point in a path on the kagome lattice. Surrounding a b step the colouring alternates 
(e.g. ABA), as in a kagome hexagon, whereas surrounding an s step it cycles (e.g. ABC), as 
on a straight line of kagome. Let {ji,j2, •••} be the steps in the loop tagged s; it can be shown 
the colours are consistent around the loop if and only if J2mi~^)"^'^'^"^ = O(mod 6). H 

Colourability does not depend on just the loop shape, but also on the orientations of the 
depleted dimers touching each site of the loop. Consider 6-loops (the shortest that we may 
have) in one layer (like those shown in Fig. [1]): each vertex of the hexagon can be tagged by 

(if the depleted dimer runs vertical from it) or 1 (if the depleted dimer sticks outwards from 
the hexagon center). Of the 2^ possible 6-loops, 45/64 are not colourable: the seven classes 
(111000), (111010), (110110), (111110), (110000), (110100), and (111100). The six colourable 
classes account for 19/64 of the 6-loops. 

Lawler's state lacks "weathervane" modes which were sometimes considered as the key feature on the kagome. 
* Note that given step i there is always (at least) one non-depleted direction available at both steps i + 1 and 

1 — 1; if the path enters on one of these but does not leave on the other, it is "bent", otherwise it is "straight". 

® Or you may just start colouring (e.g. ABAC ) in accordance with the rules mentioned, and see whether the 
loop ends on the same colour it started with. 



3. Classical cooperative paramagnet 

I first review known properties of random ensembles with divergence-like constraints. Say some 
bonds of a bipartite ^-coordinated lattice are coloured (with the same colour) such that z' bonds 
are coloured at every vertex. Construct a divergence-free vector field by letting every bond have 
unit flux from the even to the odd endpoint if coloured, or flux —z'/{z — z') if uncoloured; the 
local volume average defines a coarse-grained "polarization" P(r) with V • P = [231 [2^ I25j . 
In a random ensemble of such colourings, the entropy density behaves as 

(t(P) = (To — -a|P|^ + higher order (1) 

Combined with the divergence constraint, ([T|) implies long-range, pseudo-dipolar correlations: 

(P.(O)n(r)) oc (2) 
a|r|" 

in d dimensions. Correlations of the physical degrees of freedom are generally proportional to 
([2]). Dimer models |23j are the case z' = 1; ice models are the case z = 4, = 2. A z-colouring 
produces z — 1 flavours of P(r) (one for each colour, minus one linear dependency) hence P(r) 
in that case is a d x (z — 1) tensor. The ground states of classical spins on a bisimplex lattice - 
our original problem - satisfy the constraint X^iea ~ the polarization in that case is a 

d X m tensor where m is the number of spin components, and here too we expect ([2]) describes 
the correlations [81IM1I25]. 

The problem at hand is a depleted pyrochlore lattice {z = 3) - periodic or random - either 
with classical m = 3 component spins, or else a 3-colouring: thus, P(r) has spin or flavour 
indices. H But I'll work out, instead, the simpler case of a dimer covering on the depleted lattice; 
similar results are expected for the real problems. 

Let the dimer pattern's polarization field be P(r) on the depleted lattice {z = 3), or E(r) as 
a covering of the original {z = 4) lattice; and let D(r) be the polarization field of the {z = 4) 
dimers defining the depletion pattern. It is trivial to check that 

P = ^(E + iD). (3) 

Now, the random depleted lattice is a "quenched" model in which first D is fixed according to 
a single-component ensemble, i.e. Eq. ([T]) with a|Pp — > a^|Dp. But consider, just for now, a 
random ensemble with two (equivalent) colours of non-overlapping dimers described by D and 
E. By cubic lattice symmetry as well as D E symmetry, the joint entropy density must be 
isotropic to quadratic order: 



HE ^DE 

° 2 



-a°^(|D|2 + |E|2 + 2AD-E). (4) 



where < A < 1 is expected. This also must be the entropy density for E conditioned on a 
frozen (not untypical) D(r). Using the change of variables Q, the entropy density is ([T]) with 
P ^ P = P -|- AD a = ^oP^ and A = |('^ ~ 5)- Thus, P is a linear combination of two 
independent fields P(r) and D(r), each having correlations like ([2]), with respective prefactors 
1/a and On the other hand, a periodic depleted lattice has a uniform D(r), offsetting 

P(r) by a constant. Note that, unless a depleted lattice has cubic symmetry or small D, the 
coefficients of PaPb terms become anisotropic, (deriving from 0{DDEE) cross terms in (|3])). 



The flavour case is the large-S quantum model at temperatures low enough for coplanarity, but high enough 
to ignore the small energies selecting specific spin patterns [24] as in Sec. 12.21 



4. Conclusion 

I have shown there exist many depleted lattices (Sec. [T]) and given a picture both of the ordered 
ground state (in Sec. E]) and the cooperative paramagnet state (Sec. [3]). In the latter state, the 
disconnected correlation function of P(r) has the pseudodipolar form ([2]). Most of the results 
apply both to regular and random depleted lattices. Random depletion is an inviting way to 
include just enough constraints to maintain many features of uniform lattices, e.g. satisfying 
every triangle and allowing coplanar ground states, which are violated in rare places in the case 
of unconstrained site dilution 
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